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Introduction 



The objectives of this study are to develop a numerical 
method for estimating tidal constituents in the open ocean using 
fundamental theoretical concepts, to adapt the method to com- 
puter use, and to apply it to the analysis of the tidal regime 
in the Gvdf of Mexico. 

The ability to predict tidal behavior in the open sea is 
usefxd in oceanographic research and coastal engineering. It 
permits removal of tidal effects from deep-sea observations 
and prediction of the effect upon the tides of alterations of 
the coast line. 

Present tidal knowledge is largely limited to the results 
of harmonic analysis at selected reference stations and anpirical 
lag and amplitude corrections for closely adjacent stations. This 
technique is, of course, limited to the coast, but the work of 
Dietrich , subjectively fitting a two-constituent (M2 S' Kl) 
model (see below) to island and coastal date between coastlines 
extends the resiilts into the open ocean. The accuracy of 
Dietrich's work is considered good in dense data regions, such 
as the North Atlantic j but in sparse data areas no general 
agreement has been reached concerning tidal regimes. The dis- 
tortion of tides in shallow water cannot be accurately assessed. 
Th\is the restdts of the method are, to some degree, in error. 

A more satisfying approach is to apply the fundamental 
hydrodynamic equations and calculated luni-solar forcing function 
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to as realistic an ocean basin as possible. 

Two methods of attack are the integration of the basic 
eqiiations and the approximation of the equations by finite dif- 
ferences, leading to a solution by numerical methods. In either 
case, a boundary value problem resvilts In which Increased real- 
ism is accompanied by increased mathematical difficvdty, where- 
as oversimplification can render the problem physically meaning- 
less. In this paper discussion will be confined to numerical 
approaches. Hansen used finite differences In his treatment 
of the M2 tidal constituent in the North Sea with excellent 
resvilts, both in tidal heights and currents, as shown in Figure 
1. Friction was included and the beta-plane approximation was 
used. Rossiter's approach Involves the use of spherical 
coordinates and frictionless flow In a meridionally bovmded 
ocean of constant depth. Coriolis effects necessitate spec- 
ial care near the poles and the parameters were modified to 
insvire convergence in these areas. In another application lat- 
itudinal boundaries were imposed to approximate the Atlantic 
Ocean. Relaxation techniques were employed successfully by 
both Hansen and Rossiter. Their results and Rossiter’s recom- 
mendation that the method be applied to irregular shaped basins 
were Important motivations for the present study. 
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Tides in the Giilf of Mexico 



The tides in the Gvdf of Mexico are modest in range, but 
diverse in character. Along the Gulf Coast the tides change 
abruptly and assume many forms. Figtire 2 illustrates the 
distribution of these regimes. Along the entire Atlantic coast 
of the United States tides vary widely in range, but there is 
no variation in semi°>diurnal character. 

In the light of this constancy the abrupt changes in the 
Gulf are very dramatic. One example of considerable variation 
of the tidal regime is in the northeast portion of the Gulf . 

At Cedar Keys, a mixed, predominantly semi-diurnal regime 
exists. The tidal curves vary markedly in both directions along 
the coast. At Pensacola, 200 miles west of Cedar Keys, the 
diurnal influences dominate. Similar changes are encountered 
in other areas along the Gulf Coast. A comprehensive descrip- 
tion of the tides along the Gulf Coast and the criteria for 
regime classification have been given by Marmer (j^j . Many 
explanations have been advanced for these shifts in the basic 
form of the tides. Harris postulated in 1897 that the semi- 
diurnal tides were the extension of the Atlantic regime, while 
the diurnal tides resulted from the natural period of the Gxilf 
basin. 

Grace compiled the various theories and applied an 
inventive analytical technique in 1932. He divided the Gulf into 
five sections in each of which the depth was treated as a 
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constant. The linearized* inviscid hydrodynamic equations of 
motion were solved within each section. Continuity was satis- 
fied at the interior boundaries between sections and the bound- 
ary conditions of no normal fltix were met at the coasts. The 
studies resulted in a unique set of corange and co tidal lines. 
Based on this work Grace put forth a theory of tidal dynamics 
in the Gulfs the semldivirnal tidal wave enters the Straits of 
Florida, moves around the Gulf coast in a contra solem amphi- 
dromy and exits through the Yucatan Strait six hoxirs later; 
the diurnal tidal wave enters the Straits of Florida and prop- 
agates in such a way as to arrive at most points along the 
coast nearly simvdtaneously . 

Several topographical features of the Gulf may well have 
discernable effects on the tides. The continental slope off 
the west coast of Florida is very steep and could act as a 
reflective surface to an incoming tidal wave. The Intrusion 
of the Desoto Canyon toward the northern coast and the 
Mexican Basin are both relative deeps with unknown effects on 
the natural period of the Gulf. The Florida Keys and the 
broad shelf off the west coast of Florida could mask Influences 
of tidal waves entering the Straits of Florida. 
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3 . The Theoretical Development 

Introduction 

To describe the tides utilizing the basic hydrodynamic 
equations one must resort to describing the separate tidal 
constituents. This limits each solution to a particular known 
period and provides for a complete solution through a summation 
of constituents. Except in specialized cases this summation 
limits the problem to linearized equations of motion. The work 
described below is directed towards a solution for the principal 
lunar semidiurnal constituent, M2. Coastal observations indicate 
that knowledge of this constituent shoxild permit explanation of 
these observations and description of regimes in the open sea. 
Harmonic analyses presented by Marmer indicate that M2 and 
Kl, the lunisolar diurnal constituent, are the dominant constit- 
uents. The variability of the amplitude of Kl is from .59 
feet to ,09 feet, ignoring the results of short-term observa- 
tions. The amplitude of M2 varies from zero to 1.07 feet by 
the same criterion. Thus the spacial distribution of the M2 
constituent should be indicative of the distribution of the forms 
of the tides. 

The Equation of Motion 

From the eqxxations of motion and continuity we may derive 
a second-order elliptic partial differential equation relating tidal 
height, depth, and their derivatives. The analysis is performed 
on a beta-plane, Advective and frictional terms in the equation 
of motion are neglected. The hydrostatic approximation is made. 
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The beta-plane assumption restricts the size of the 
basin in which the res\ilts will apply, but is considered accep- 
table in the Gulf . 

Friction is present and does influence tidal distributions, 
especially in shallow ■water. Frictional effects are, however, 
believed to be highly non-linear, thus further complicating the 
equations. (The work of Dronkers(|^ in linearizing this term 
for shallow water is an example of these complications.) Fur- 
ther justification for neglecting friction Is that only the per- 
imeter is in shallow water while most of the Interior Is in 
relatively deep water. 

With the stated assumptions the dynamic equations are; 



where, 

u is the velocity component In the x direction (positive 
East) , 

V is the velocity component In the y direction ( positive 
North) , 

F Is the Coriolis parameter, 
g is the acceleration of gravity, 

^ is the tidal height above vmdisturbed sea level (which 
is mean low water) and is positive upwards, and 
^ Is the forcing function (see below). 

The hydrostatic approximation implies that u and v are 




( 1 ) 



Independent of depth. 



( 2 ) 



The eqxiation of continuity iss 












•where h is the basin depth. 

Since the frequency of the desired harmonic constituent 
is known, sinusoidal forms for the dependent variables are 



assumed: 

u- < 3 ) 

H Cos O^t 

where the notation, Re ^ J indicates the real part of the 
quantity in brackets, 

a is the angular frequency of the constituent tinder exam- 



ination, and 

H is the amplitude of the forcing function. 

Solving the equations of motion for u and v and substitut- 
ing the results into the continuity equation yields one partial 
differential equation for the tidal height. Writing 
(where, for instance,^ is the component in phase wi-th the 
forcing function) and separating the differential equation into 



real and imaginary parts gives 

h(K-n) _ F. ^ 



p. ah. a A 4. f :=■ o 

■F O' ^ 3hT J# - ^ 



O' a x 






and 



. F ah 0-P jj = 0 

<r 57 



X. 



a 



( 4 ) 
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•where 



il ^ 

Rossiter used spherical coordinates In treating basins 
covering large areas of the ear-th. Thus his eqiiations differ 
in detail at this point although the form is very similar. 

Hansen's work In the North Sea Included friction. The 
added terms in the eqviations of motion had the form, 



F^-- R and ’R'J (5) 

where R Is a suitably chosen proportionality factor. The 
excellence of his residts is a strong argument for this some- 
what ad hoc hypothesis. In the situation involving considerably 
greater bottom relief this form is inadequate due to the var- 
iability of the friction effect with depth. In deep water 
friction will not Influence tidal distributions as much as bottom 
topography which enters the equations through continuity con- 
siderations . 



Boundary Conditions 

The boundary conditions determine the velocities normal 
to the boundaries. At solid boundaries these velocities must 
be zero. At the passages Into the Gulf values must be assign- 
ed to the associated tidal currents. Following the treatment of 
tidal heights, we write 

U = UKI + L UK2. 

V - VKl I VK2-. 



( 6 ) 
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Substituting these relations into the solutions for u and v fr*om 
the equations of motion and separating into real and imaginary 
parts gives 



Observations of tidal currents to determine their con- 
stituents are not frequently made., For accuracy continuous 
observations are required over long periods of time. Tidal 
cxxrrent measurements made by Harris In 1887 and supported 
by Grace's work were used in the Yucatan and Florida Straits 
in this study. The accuracy of these calcvilations is question- 
able but the order of magnitude is correct. By allowing flow 
through these channels the effects of a co-oscillating tide are 
included in the solution. 



The forcing fiinction Is basically the gravitational effect 
of the ideal body associated with the constituent to be examin- 
ed, It is obtained from equilibrium tidal theory (for example 




The Forcing Function 




17 



for the Ivinar semidiurnal constituents, and 



,0d^ H^Y^ 2 f c^s 'I <K>S 

H - CL (gj \_ 4 




for the l\inar diurnal constituents. 

Here, 

a Is the radius of the Earth 

d is the distance from the Earth to the Moon 
Mm is the mass of the Moon 
Me is the mass of the Earth 



is the longitude relative to the reference longitude. 

For computations the central meridian of the Gulf, 90°W, 
will be used as a reference. 

Similar expressions apply for solar diurnal and semidiurnal 
constituents and also for those with periods of up to two weeks. 
In this work these constituents have relatively small amplitudes 
and so are not treated. 

Over the area of the Gulf the lunar semidiurnal forcing 
function varies by about eighteen percent. 




is the latitude of the point under consideration 




is the declination of the ideal body 
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4. Numerical Techniques 

Finite Difference Approximations 
To put equations (4) into a form more suitable for num- 
erical treatment the derivatives are approximated by central 
differences, 



Z6 



, fIfV' 






and 



I ax Jo 






( 10 ) 



( 11 ) 



where2 is a dummy variable and the subscripts refer to loca- 
tion in the grid arrangement as shown below, _ 

Yvl 

f[ 



y-i 



x+i 






Central differences utilize values sampled at an interval of two 
grid distances. This has a smoothing effect on the data. 

Because the bathymetry is normally given to the nearest fathom 
smoothing is beneficial in decreasing errors in depth gradient 
due to sotmding inaccuracies. The tidal heights are much smaller 
than the depths. Because of this, solutions are very sensitive to 
sounding errors. For this reason central differences are used 
in all finite difference approximations. 

The use of central differences at a point on the boiind- 
ary introduces values for tidal heights outside the boundary into 
the equations. 

For an interior point, equations (4) with finite differences 
become 
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-t- ('4ho-Hx+, + Wy-;) 

+ (ST-I^)^^y (a Uo+ Y-'^ ■*■ W-1 (^4 Uo "Wy+i + ^v-'') 

+ ^ [|: (^Y« ■'^Y->)3 +(^)y», 

- (^)o ^ 



and 



4- (^4Uo-Kx^ i 

(-^)y+\ ("41no Wi^-h ” ^4-1^ ■^(^)y.\ (4ho’ hY+v"^ Iav(m'^ 

+ ^"H)x4 (" /T-w^y-^-i 

f ^ (-^Jo ^o. 



Consideration was given to variable grid distancing to permit 
detailed analysis In areas of large depth gradients. However, 
the stability of the computation about a point with different 
distances to each of the adjacent points Is questionable, and 
the finite difference form of the equations is quite complex. 
Appendix B contains the finite difference equations for a 
variable grid distance. The density of bathymetric data In 
the Gulf Is such that a grid finer than that used below woxild 



20 



require considerable interpolation and thus no improvement in 
the solution would appear. 

Rossiter used a finite difference equation involving der- 
ivatives at the boundary and three interior points to produce 
the tidal heights outside the boundary. Although it is valid, 
Rossiter's procedxxre has the disadvantage of Introducing para- 
meters associated with non-adjacent points into the equations. 
This complicates the system of equations In a relaxation tech- 
nique by Increasing the number of points involved in each equa- 
tion by two. It also disrupts the iterative form designed for 
computer use, which will be discussed later. 

At a boiindary point, the basic equations (12), derived 
from continuity considerations and the applicable boundary con- 
ditions are both in effect. These boundary conditions manifest 
themselves in the forms of prescribed velocities normal to the 
boundaries. 

The velocity equations (7)in finite difference form are 



_ - < 3 ' 

VK 1 ■ 2.(0*- f'7 

\Ikz-- 









UKZ- 






(13) 
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Thus, for a given normal velocity, the tidal height for the 
point outside the boundary may be solved for in terms of three 
Interior tidal heights. Substituting these solutions into the 
equation for an Interior point now incorporates the dynamic 
equations, the continuity equation and the boundary condition 
Into a set of two eqxiations at the point, as well as eliminat- 
ing values outside the boundary from the calculations. 

The substitutions and resultant equations will be unique 
for each bovindary or corner. As an example, to solve for a 
point on the south boundary with a known velocity, equations 
(13) are used in the form. 



and substituted into the Interior equations ( 12 ) to eliminate 

values Involving the exterior point Y- 1 . 

This results in 



+ rjT-n)x», + 




(14) 






and 
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+ (j?) ^ ^)^+‘ " W/-»)(\ “ 3 



+- VK\ ^0"“ 

- VKiP*f|E^‘'(h>,-h,-0] --0. 



To accommodate the boxmdary conditions involved at a corner 
values of u and v must enter the basic equations (12). This 
is accomplished by simtiltaneous solution of equations ( 13 ) to 
eliminate terms involving the two heights at the same point 
from the combined equations. The results are 






/ ^\ / ^ N T26 FC<rl^ 1 1 

(JiJxtl " ( ^zJx-i + VKZ Lat i ( I - p^^») J 

( 16 ) 



Combinations of these equations will meet any corner boundary 
condition. For example, substituting to eliminate values at the 
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point X+1 in the south boundary equations will give equations 



for the southeast corner in the form, 



u 



+ (JT-^)y 4| ^8(^o] + VK\ ^ U 

- UK2. ■^)r-^. + hx4rbA-,)]] 



- VKZ. 



[I- 



■h 



4^A CQ--F'-') 



2 6 F 






- 0 



)(4-ho- 



(17) 



and 



+ (J1)y+i + VK2. C^C'^Ci-p‘Ar»'^ ^ D 

■*• U K I )( 



■^■3- 






)(4loo-h 



lO )(4-t ■" 

, ^ ho \1 

i>^+i ho\-» ■*■ cT^( V - 



VKl[X' 

-o. 



The forms of the equations for each boimdary and corner 
point are presented in Appendix A, 
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Methods of Solution 



When eqxiations (12) and their equivalents incorporating 
bovindary conditions (see Appendix A) are applied at each point 
of a grid consisting of n points describing the body of water 
studied, a system of 2n equations in 2n unknown tidal heights 
results. Both methods of solution undertaken here are pred- 
icated upon the presence of a dominant diagonal in the system's 
coefficient matrix. That is, in the application of the equations 
to each point there is one coefficient much larger than all 
others in the equation and in each successive equation the 
largest coefficient multiplies a different ixnknown value of the 
tidal height. Both methods are designed to obtain a solution 
by altering the unknown tidal height associated with the largest 
coefficient. 

Gauss - Seidel iteration (see for example, Hildebrand CsD ) 
is the method used in the computer program. Briefly, this 
method consists of obtaining new values for the variable with 
the dominant coefficient, say , in an eqxiation of the form, 

A| X| + A 2 . X-t + • • • -f Am Xh “0 

by the formula 




This satisfies the equation with approximate values of X^,...., 
Xn by varying X^ only. This procedxire is carried out for each 
equation in the system following a first guess of the value of 
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f 



each variable. Repeated use of equation (18) usually brings 
about "convergence", in the sense that each of the equations 
is satisfied by the same set of unknowns to within prescribed 
tolerances. Gauss - Seidel iteration is not highly sensitive 
to the first guess values. 

Relaxation is a similar process but permits selective 
variation of the unknowns to bring about convergence. Using 
the same example equation the first guess will give 

A\ X\ Az)^ 2. + • • • AinXo “ 

where R Is the "residual", which should vanish to achieve con- 
vergence of the system of eqviations. Changes will be made 
in the unknown associated with the largest coefficient following 
the rule 

x/ = Xi ■ 

where 7^ is the "Relaxation Coefficient" and may take on values 



in the range 



0 < > < Z. 



( 20 ) 



■Whenl^ =1 the residual will be reduced to zero and the 
process is very similar to the iteration described above. When 
7'>| the residual tends to alternate sign on successive iterations 
and decrease in magnitude. This process usually speeds up 
convergence and Is called over relaxation. WhenT'^lthe residual 
tends to be reduced in raagidtude and keep the same sign thus 
approaching zero from one side rather than oscillating around 
zero. This is called underrelaxation. 
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5o Hand Calculation of the M2 Tidal Constituent in the Gulf 

To gain an understanding of the behavior of a tidal solu- 
tion in the Gulf and to supplement the computer •work a hand 
solution was undertaken for a relatively coarse grido The grid 
arrangement consists of 16 points with a grid distance of 180 
nautical miles Figure 3 shows the arrangement. This grid 
represents the Gulf fairly well, but has the disadvantage that, 
in areas of large depth gradient, the terms in the equations 
may not occur with the relative magnitudes desired for a rel- 
axation solution. 

Applying equations (12) and appropriate boundary conditions 
to each of the l6 grid points produces a system of 32 equations 
for 32 \inknown tidal heights. Of these, eight do not contain 
the dominant coefficients desired for relaxation. Six of these 
eight are for interior points wi'th eight terms in each equation. 
This means that successf\il reduction of the residuals of these 
equations, especially in the later stages of the solution, may 
depend upon systematic adjustment of the number of tidal heights 
to avoid deleterious effects on the residuals at adjacent points. 

In every iteration values of the tidal heights used were 
the results of the last calculations made at the point. Some 
systems (e.g., the method of Richardson as described by 
Hildebrand [i b!d^ ) are cyclic in nature using only calculations 
from pre-vious iterations. The approach in this work is not 
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cyclic, but is based upon maximum reduction of the largest 
residuals consistent -with the long-range effect upon adjacent 
points. 

Initial passes were made based on a first guess field de- 
rived from Grace’s results and observational data reported by 
Marmer. Over-relaxation techniques were utilized with tak- 
ing values from 1,5 to 1.2. The chief drawback of this tech- 
nique is that It can produce oscillations In the values of the 
tidal heights and thus fail to converge. At the first sign of 
an undamped oscillation,*?^ was lowered to 1. On subsequent 
Iterations the relaxation factor was varied to meet the needs 
of each equation. To determine corrections late In the cal— 
ciilatlons systems of simultaneous equations were solved for 
the appropriate values of *1^ . At this stage the equations 
were very sensitive, and large corrections frequently had del- 
eterious effects on adjacent points. To most efficiently deal 
with these situations tabulations of present resldvials and the 
effects of proposed changes In tidal heights upon them were 
necessary. These were most helpfvd In determining the Inter- 
action of changes In the tidal heights. After twenty Iterations 
component tidal heights (^&nd/g) at all Interior and solid bound- 
ary points were defined to an accuracy of plus or minus three 
centimeters. Those at the Yucatan Passage and Straits of 
Florida were extremely sensitive to values of Input cixrrent, 
fluctuating approximately ten centimeters with changes of only 
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one centimeter per second in the associated velocity component. 
The effects of these points on other equations is not large. 
Since the current velocities utilized are not considered suffic- 
iently accurate to warrant it, these points were not relaxed 
to the degree applied to the other points in the grid. 

The resxilts of these calculations are contained in Appendix 
D, In general, this solution was very similar to that accom- 
plished using the computer on the same 16 point grid. The 
maximum difference in magnitude of tidal height was 9_»0S cms. 
The average difference was 3.23 cms. This computer solution 
is discussed in more detail below. 



29 



6. Solution by Digital Computer 

Adaptation of the finite difference formulae (12) to the 
computer is a difficult task. For efficient iteration a proced- 
ure is necessary for expressing all equations in a common for- 
mat. This procediire permits the use of indexing operations 
in the identification and location of equations. This utilization 
of indexing is accomplished by the inclusion of identification of 
adjacent points and their relative positions in the input data. 
Throixgh this, and a coded method for points outside the bound- 
ary the program identifies interior and specific boundary loca- 
tions and simiiltaneously applies the proper equations to them, 
Fortran 60 program langiiage is used. Based on the needs of 
this particiilar application, convenient computer capacity, and 
the concept of versatility, the maximum number of points in 
the grid is 1200. These points are numbered serially, and 
points outside the boimdary are identified by numbers greater 
than 1200. Since the reqviirement for ficticious tidal heights 
outside the boundary has been eliminated (see above), the only 
data required outside the boundary are tidal ciirrent velocities 
and basin depths. 

The presence of an identification number greater than 
1200 then specifies a particiilar boiindary and two such numbers 
specify a particular corner. The cxirrent velocities included as 
data for points outside the boundaries replace tidal heights for 
those points in the equations at the boundary point. This 
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req\iires alteration of the coefficients in the equations to fit 
the location of the point at which the equations are being ap- 
plied. This is accomplished from the general forms, 



Qa, 


- b, 




^ Jz) x*i 

jx (21) 




+ b3 C 


0. 4 ('^-^) Y+\ 
+ cj/ CZ 

and 


- b4 ( ) Y* 1 




f bi 


^CX-i ( fl) 


- bi ( K-i 




- bs f 


^ 


+ b 4 ( "S'' "M) y-v3 / C 2. . 



Thtis the coefficients depend upon the location of the point. 
Appendix A presents these coefficients in detail and the com- 
puter program is given in Appendix C. 

The system of equations derived from the grid utilized 
in the hand calculations was then solved using the computer. 
Actual depths and currents were employed. This system con- 
verged in 37 iterations. 

In a study of this nature results are best displayed by 
the analysis of the solution to obtain cotidal and corange lines. 
Gotidal lines are lines along which the high tide occurs simul- 
taneously. Corange lines are lines along which the tidal ampli- 
tudes are equal. As set forth above these resxilts apply to the 
M2 constituent, and the phases are relative to the passage of 
the moon through the 90th west meridian. The lack of detail 
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In this solution causes the corange and cotidal analysis to be 
somewhat subjective. The analysis of the 16 point solution is 
illustrated in Figure 4 which is based upon coastal data as well 
as the computer solution. It is evident that most of the 
distortion of the tides occurs in the coastal regions. This 
does not mean that the solution is in error. The lack of 
detail does not permit identification of a small-scale amphidromy 
such as that proposed by Grace or seen in the solution discussed 
below. The path of the progressive tidal wave indicated by the 
cotidal lines Is logical and amplitude distributions agree with 
observed coastal data. 

To obtain more detail in the solution a finer grid was 
desired. Using a grid distance of thirty nautical miles a grid 
was then constructed in an irregular shape closely approximating 
the Gulf. Real depths were used, taken from U. S. Coast 
and Geodetic Survey Chart 1007. The grid consisted of 504 
points. Current data mentioned above were used in the Straits 
of Florida and Yucatan Strait. The resiJlting system of equa- 
tions failed to converge after 2500 Iterations. In the initial 
iterations the tidal heights increased slowly in magnitude. 

After 25 Iterations the solutions were unrealistic and they 
Increased more rapidly in magnitude. Computational Instability 
seemed to be present and further efforts were directed at a 
more tractable situation. 

Using the same grid a constant depth of 1000 fathoms 
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was introduced. The system of equations converged in 378 
iterations with a relative convergence criterion of 0.1. That 
is, each tidal height within the grid must change by less than 
10^ of that of the previous iteration before the iterative 
procedure is stopped. The maximum change of tidal height 
over the last three iterations was .12 centimeters. These 
results are tabulated in Appendix D. Figure 5 shows the 
analysis of these results in terms of corange and cotidal lines. 
The path of the M2 constituent tidal wave crest through the 
Gulf is described by the cotidal lines labeled in hours with 
zero being the time of the passage of the moon through the 
90th west meridian. 

The detail obtained in a 50k point grid considerably lessens 
subjectivity in these cotidal and corange analyses. Significant 
aspects of these analyses are: the well-developed amphidromy 
in the Northeastern corner of the Gulf, the presence of the 
greater amplitudes arotand the coast, the presence of a rather 
flat field of low amplitudes in the vicinity of the amphidromic 
point, and the relatively high amplitudes in areas of known semi- 
diurnal dominance. The path of the tidal wave crest through 
the Gulf is very well defined. 

Both the 16- and 504- point grid computer solutions 
produce corange and cotidal lines consistent with coastal data. 
Because one is based upon constant depth and the other upon 
actual depths the analyses differ considerably in some areas; 
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however, the regime shifts mentioned in Section 3 (say, 
semidiurnal to diurnal) are clearly indicated in both Figures 
2+ and 5« 
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7. Conclusions and Acknowledgements 

The three solutions obtained indicate that the numerical 
methods previously discussed are effective in predicting tides 
at sea. Based on the relatively quick convergence with actual 
depths of the 16 point system which has a grid distance of 
180 nautical miles, the convergence with constant depth of 
the 504 point system which has a grid distance of 30 miles, 
and the failure to converge with actixal depths of the 504 
point ' system, it is reasonable to assume that there exists an 
intermediate grid distance such that convergence may be ob- 
tained with maximum detail in the solution. Another approach 
is to obtain the maximum detail within the field by mathemati- 
cally smoothing the actiial depth data to such a degree that the 
504 point system would converge. The danger here is that 
excessive smoothing could distort the data thus giving less 
realism in the solution. 

The cotidal and corange analyses of the constant and 
variable— depth models indicate that variations in depth are not 
solely responsible for the tidal regime shifts briefly described 
in Section Is (see Figure 4). Additional studies may deliniate 
the relative importances of spacial forcing function variation, 
input current, basin shape, and basin depths upon tidal ampli- 
tude distributions and perhaps uncover other important variables. 

By removing certain topographical featiires of the Gulf floor 
from the data it would be possible to assess their effect upon 
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tidal distributions. Similarly variations In input currents and 
the forcing function field could be examined for their effects. 
By adjusting Input currents only, under a condition of maximum 
realism and detail in all other data, to obtain the constituent 
distribution most compatible with coastal observations one would 
be able to evaluate the accuracy of present constituent tidal 
cvirrent data. 

The only geographical limitation In this technique Is the 
size of the basin due to the beta-plane assumption. Thus all 
other seas comparable in size to the Gulf of Mexico are well 
suited for an investigation of this sort. 

This study and the writer have both benefitted greatly 
from the Interest, guidance and considerable abilities of mem- 
bers of the facvilty and staff of the United States Naval 
Postgraduate School. Principal among these are: Associate 
Professor T. Green, III for his encouragement, help, con- 
siderable effort and constructive guidance; Professor G. «J. 
Haltiner for his counsel and personal example of scientific 
excellence; Associate Professor *J, B. 'VN^ckham for his Incisive 
and constructive criticism; R. D. Brunell for his personal 
Interest, efforts In the programming and technical assistance; 
also Dr. G, B. Austin and Dr. G. Dowling of the U. S. Navy 
Mine Defense Laboratory, Panama City, Florida for their 
Interest and assistance In selection of this topic; and Mrs. 
«J^alene D. Gainer for her patience, encouragement, understand- 
ing and recently acquired proficiency at the desk calculator. 
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Figur© 3. - GRID ^ARRANGEMENT FOR HaND 
CALC UL ATIO N S 

This consists of l6 point grid with ISO nautical 
mil© grid distance. Input currents war© u at point 14 
and V at point 13 . 
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Figure 4. - RESULTS OF VARIABLE DEPTH 
COMPUTER SOLUTION: 



Corange lines (dashed) and cotidal lines (solid). 
Corange lines are labeled in centimeters. Cotidal lines 
are labeled in degrees denoting the phase of the tide 
relative to meridional passage of the moon through 
90OW. 
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Corange lines (dotted) ‘.abelod 
Cotidal lines (solid) labeled in hours 
time of the moon's passage fnrough 
meridian. 



in centimeters, 
with sero being 
the 90th West 



the 



APPENDIX A 



FINITE DIFFERENCE EQUATIONS FOR INTERIOR 
AND ADD BOUNDARY POINTS ADAPTED FOR 
COMPUTER SODUTION 

Equation (12) represents the finite difference form for 
an interior point. We have shown that at each boxmdary un- 
ique modifications of this equation are necessary to represent 
continiiity and the boundary conditions. To avoid the necess- 
ity of separately establishing each equation within the program 
and to create a compact iteration scheme the following system 
was devised. 

Examination of all the eqioations led to a general form, 
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are divided by C2 
and 
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These equations are set up to carry out Gauss-Seidel 
iteration, but modification for other solution methods could 
be made at this point. 



Now by indexing within the program the constants may 
be called to fit any of the previously derived eq\xatIons satis- 
fying both dynamic considerations and boundary conditions. For 
interior points the constants are: 
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APPENDIX B 



DIFFERENCING EQUATIONS FOR VARIABLE 
GRID DISTANCE 



Given a relative grid arrangement of the sort, 



Y-^-\ 




it; is desirable to use only contiguous points in the differencing 
scheme because the use of other points would require a field 
of values and points outside the real grid. This in turn would 
negate the use of central differences in evaluating the second 
derivative. 

The basic forward differencing formulae were applied to 



7 _ (AX-) >:4-v) 4- - 7o 

VC- 

, (aY-) - 7o[(AYi) 4 -(aY-) ~] 

(aY-)^ AYj. 



where H is a dummy variable. This expression is qviite involved 
and the stability of calcvilation about such a point is highly 
questionable. 
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PROGRAM TIDES 

dimension N( 1200*A) » con (1200 *10) »FF( 1200) *ZL (1200) »HF( 1200) . 
1WL( 1200) *0EPTH( 1401) *2 ( 1401 ) *W( 1401 ) 

COMMON N » CON ♦ FF ♦ 2L tHF ♦ WL » DEPTH *.2 » W 
EQUIVALENCE ( FF »2L) ♦ ( DEPTH*WL ) 



APPENDIX r 

THE COMPUTER PROQRAM 
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C0N( I ,6)=C0N( I ,5) 
C0N( I ,7)=A2*CXC 
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110 ACTAN=-PI+ATANF(A/B) 230 

200 RETURN 240 

END 

END 250 



APPENDIX D 



THE 


NUMERICAL RESULTS 
CONSTANT DEPTH 


OF THE 504 POINT 
SOLUTION 


POINT 


MAGNITUDE (cms) PHASE ( 


1. 


1.3887000 


- 124 . 1.245 18: 


2. 


1.3756521 


-I38.8O64I 


3. 


1.3217907 


-131.90854 


4. 


1.2539393 


-120.89377 


5. 


1.2029531 


-105.43202 


6 . 


1.2059750 


-86.843942 


7. 


1.2334819 


-73.334869 


8. 


1.7517582 


-172.50711 


9. 


1.7194041 


-172. 30618 


10. 


1. 7070080 


-171.97713 


11. 


1.6958950 


-171.56044 


12. 


1.3027698 


-140.26361 


13. 


1.2262120 


-137.87719 


14. 


1.1332772; 


-130.61788 


15. 


1.0404123 


-118.38878 


16. 


.98849004 


-99.810982 


17. 


I.07I8649 


-.75.709188 


00 


1.4355376 


-53.916750 


19. 


2.4490565 


-39.475631. 


20. 


4.1602358 


-23.200350 


21. 


7.4439120 


-4.1130445 
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22. 


5.4827457 


-2.1818875 


23. 


1.6776942 


-172.36824 


24. 


1.6361282 


-172.15673 


25. 


1.5939419 


-171.76625 


26. 


1.5654488 


-171.39348 


27. 


1.5564011 


-170.96102 


28, 


1.5592573 


-170.38945 


29. 


1.5838993 


-169.61342 


30. 


1.5963884 


-167.86595 


31. 


1.5907292 


-165.26372 


32, 


1.5549091 


-161.59242 


33. 


1.4619722 


-156.27691 


34. 


1.2553052 


-148.15227 


35. 


1.1279148 


-142.15830 


36, 


1.0039504 


-134.12128 


37. 


.88153194 


-131.44064 


38. 


.79735342 


-100.55519 


39. 


.86428115 


-71.147404 


40. 


1.2173230 


-45.294358 


41. 


1.9446390 


-28.377334 


42. 


2.9908028; 


-15.251722 


43. 


4.1469681 


-3.9740857 


44. 


3.7697796 


1.7862483 


45. 


2.9949950 


9.8396774 


46. 


1,5900268 


-171.85187 
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47 



< 

49. 

50. 

51. 

52. 

53. 

54. 



55. 

56. 

57. 

58. 

59. 

60 . 



61 . 



62 . 

63. 

64. 

65. 

66 . 

67. 

68 . 



69. 



70. 

71. 



1.5265542 

1.4871450; 

I.48I7022; 

1.4258745; 

1.4202776: 

1.4242894 

1.4326960 

1.4367689 

1.4264928. 

1.3910475. 

1.3153951 

1.1943604 

1.0752430' 

.94189596 

.79343153 

.64595142 

.59415964 

.80011771 

1.2755715 

1.8856044 

2.4422500 

2.5285018 

2.5386045 

1.5837556 

1.5189584 



-171. 55308 
- 171.29280 
-171.10465 
-170.47495 
- 169.98678 
-169.33126 
-168.36735 
-166.84673 
-164. 70118 
- 161.82801 
-158.02929 
-153.24541 
-148.40719 
-142.03557 
-132.02672 
-113.55350 
-79.717909 
-44.336330 
- 22.528501 
-9.0438157 

1.3021725 

10.436640 

20.977383 

-171.35588 

-171.12447 
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72. 


1.4459243 


-170.90635 


73. 


1.3873289 


-170.62673 


74. 


1.3523482 


-170.37487 


75. 


1.3203357 


-170.00174 


76. 


1.3014904 


-169.57599 


77. 


1.3007594 


-169.12103 


78. 


1.3081634 


-168.52707 


79. 


1.3162539 


-167.68322 


80. 


1.3231527 


-166.52922 


81. 


1.3207365 


-165.02452 


82. 


1.3008327 


-163.15918 


83. 


1.2537422 


-160.91726 


84. 


1.1791414 


-158.41437 


85. 


1.0871412 


-155.84092 


86. 


.9684444 


-152.64950 


87. 


.82084260 


-148.05767 


88. 


.63568331 


“139 . 73622 


89. 


.42310620 


-119.25154 


90. 


.32467912 


-66.502209 


91. 


.55898858 


-18.151299 


92. 


.93149353 


2.2351259 


93. 


1.2945355 


15.742440 


94. 


1.5414282 


29-348534 


95. 


1.8510571 


43.251027 


96. 


1.4426962 


-170.27041 
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97. 


1.3801700 


-170.11656 


98. 


1.3135581 


-169.87776 


99. 


1.2602636 


-169.60957 


100. 


1.2306716 


-169.39582 


101. 


1.2046669 


-169.06696 


102. 


1.1925103 


-168.71396 


103. 


1.1990822 


-168.38716 


104. 


1.2142370 


-167.98880 


105. 


1.2307824 


-167.44851 


106. 


1.2495502 


-166.79710 


107. 


1.2644726 


-I66.O42O8 


108. 


1.2694225 


-165.22816 


109. 


1.2565822 


-164.40579 


110. 


1.2247704 


-163.71041 


111. 


1.1738689 


-163.24897 


112. 


1.0947432 


-162.92580 


113. 


.98741067 


-162.80156 


114. 


.84121204 


-162.89192 


115. 


.64327279 


-163.50925 


116. 


.40957519 


-167.70662 


117. 


.16316727 


160.04043 


118. 


.27612851 


75.916858 


119. 


.60679389 


67.911211 


120. 


.98720807 


74.397350 


121. 


1.4865882 


81. 249956 
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1.3192293 
1.2570372 
1.1937553 
1.1450593 
1.1211817 
1.1020564 
1.0981462 
1.1145761 
1.1415354 
1.1725985 
1,2102442 
1.2500523 
1.2873929 
1.3153326 
1.3320425 
1.3344956 
1.3135021 
1,2716266 
1.2020871 
1,0981960 
.9868950 
. 88220144 
.86646759 
,98653169 
1.2799964 



- 168.98529 
- 168.94400 
- 168.73472 
- 168.49526 
- 168.34803 
- 168.10229 
- 167.88417 
- 167.77086 
- 167.67201 
- 167 . 54610 
- 167.46968 
- 167.47233 
- 167,60798 
- 167.92557 
- 168.52434 
- 169.47107 
- 170.77450 
- 172,57444 
- 175.05351 
- 178,68367 
175.47563 
165,61313 
150,92426 
136.05371 
126.90135 



147. 


1.7474118 


121.69064 


148. 


1.2092843 


-167.64318 


149. 


1.1461983 


-167.67550 


150. 


1.0849542 


-167.48544 


151. 


1.0403817 


-167.27501 


152. 


1.0224132 


-167.21450 


153. 


1.0109213 


-167.08640 


154. 


1.0166413 


-167.05710 


155. 


1.0451478 


-167.22489 


156. 


1.0872336 


-167.49308 


157. 


1.1373998 


-167.83048 


158. 


1.1995002 


-168.32647 


159. 


1.2705283 


-169.00223 


160. 


1.3470828 


-169.88820 


161. 


1.4228629 


-171.00051 


162. 


1.49577548 


-172.38357 


163. 


1.5614789 


-174.05059 


164. 


1.6100079 


-175.98955 


165. 


1.6439044 


-178.23670 


166 . 


1.6564831 


179.18471 


167. 


1.6410146 


176.18332 


168. 


1.6223832 


172.66604 


169. 


1.6032476 


168.28565 


170. 


1.6392330 


161.61001 


171. 


1.7547253 


155.14506 
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172. 


2.0045504 




150,74902 


173. 


2.3098794 




149.75370 


174. 


3.0581112 




155.39824 


175. 


1.1086753 




-166.33404 


176. 


1.0442983 




-166.35663 


177. 


.98475467 




-166.13134 


178, 


.94421575 




^165.92555 


179. 


.93236154 


c 

/f'’ 


-165.95960 


180. 


.92907261 




-165.97488 


181, 


.94545533 




—166 .17686 


182, 


.98770361 




-166,67282 


183. 


1.0474125 




-167.33990 


184. 


1.1201448 




-168.13757 


185. 


1.2111263 




-169.14457 


186 . 


1.3187259 




“I70. 35560 


187. 


1.4407064 




-171.76804 


188. 


1.5713587 




“I73.366I6 


189. 


1.7083588 




-175.15459 


190, 


1.8462118 




-177.11356 


191. 


1.9728861 




-179.21124 


192. 


2.0883663 




178,57807 


193. 


2,1825526 




176.30378 


194. 


2.2449039 




174.02875 


195. 


2.2967387 




172.02866 


196. 


2.3377767 




171.09147 
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197. 


2.4036320 




165.81038 


198. 


2 . 5284900 




162.15293 


199. 


2.7636824 




159.89052 


200. 


3.1167505 




158.85404 


201. 


3 . 5327726 




157.76749 


202. 


1.0129899 




-165.13166 


203. 


.94788364 




-164.99387 


204. 


.89049944 




-164.63719 


205. 


.85424915 




, -164.38986 


206. 


.84869031 




-164.51702 


207. 


.85389577 




-164.69366 


208. 


.88147957 




-165.15934 


209. 


.93842593 




-166. 01340 


210. 


1.0173279 




-167.08388 


211. 


1.1149818 




-168.30477 


212. 


1.2381205 




-169.73308 


213. 


1.3866221 




-171.33510 


214. 


1.5595657 




-173.08435 


215. 


1,7520227 




-174.95054 


216. 


1.9616174 




-176.92430 


217. 


2.1814858 




-178.97981 


218. 


2.3965146 




178.91438 


219. 


2.6017553 




176.79859 


220. 


2.7802197 




174.727I8 


221. 


2.9135588 


0 


172.74419 
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222 , 3o 

223= 3= 

224= 3= 

225= 3= 

226 = 3 = 

227= 3= 

228. 4= 

229. 5= 

230, Ic 



231 = 
232, 
233 = 
234. 
235 = 
236. 

237 = 

238 = 
239 = 



240 . 1 . 

241, 1, 

242 c 1 = 

243 = 1 = 

244= 

245 O 2 o 

246 e 2 , 



0161880 


170,95484 


0826239 


169=29360 


1604651 


167=06122 


2728319 


165.33070 


4899034 


164=56022 


8443713 


164.85080 


3693820 


166.24802 


2322027 


169=77120 


0182813 


-165.60600 


91611772 


-164.13332 


85362642 


-163.50245 


79976933 


-162,89356 


76826352 


-162.55455 


76900960 


-162,77557 


78258238 


-163.13041 


82125881 


-163.88891 


89298309 


-165.12379 


99155335 


-166.58995 


1152134 


-168.18580 


2724424 


-169=94772 


4649441 


-171.82287 


6935393 


“173. 77752 


9546477 


-175.78511 


2464576 


-177.84262 


56IO4I8 


-179.93537 
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247. 


2.8795196 


177.95661 


248. 


3 . 1896662 


175.87351 


249. 


3.4633445 


173.86999 


250. 


3.6690512 


171.98936 


251. 


3.8103315 


170.31615 


252. 


3.8769698 


168.82943 


253. 


3.9217889 


167.55483 


254. 


3.9790625 


166.78406 


255. 


4.1409762 


166.88418 


256 . 


4.4609625 


167.91345 


257. 


4.9953985 


169.71161 


258. 


5.8880954 


171.96057 


259. 


7.2656523 


173.31975 


260. 


8.1811000 


175.57659 


261. 


.90562088 


-163.42604 


262. 


.81687892 


-162.22822 


263. 


.76044280 


-161.48515 


264. 


.71118680 


-160.67153 


265. 


.68455013 


-160.23044 


266. 


.69116844 


-160.56392 


267. 


.71238788 


-161.12164 


268. 


.76123619 


-162.21030 


269. 


.84674989 

1 


-163.85739 


270. 


.96413752 


-165.71847 


271. 


1.1133331 


-167.65178 
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272. 


1,3048910 


-169.68291 


273. 


1.5428370 


-171.75315 


274. 


1.8304559 


-173.83722 


275. 


2,1665757 


-175.92528 


276, 


2,5513942 


-178.03303 


277, 


2,9773257 


179.83285 


278. 


3.4219972 


177.67307 


279. 


3,8637190 


175.52600 


280. 


4.2571099 


173.45832 


281. 


4.5479984 


171.53646 


282. 


4.7205917 


169.86121 


283. 


4.7548153 


168,47200 


284. 


4,7096263 


167.52009 


285. 


4.6451493 


167.17105 
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-35.224494 
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NUMERICAL RESULTS OF HAND CALCULATIONS 



POINT 


MAGNITUDE (cm) 


PHASE 


1, 


1,414 


315° 


2, 


9.86 


24° 

0 


3 « 


10,1 


354 


h « 


1.414 


45 ° 


5 « 


7 


0 


6 . 


14 


0 


7 . 


3 


0 


8 . 


11.4 


343.3° 


9 . 


6 


90° 


10 . 


3.14 


14° 


11 . 


7.07 


^ 0 

8,14 


12 . 


6.0 


0 


13 . 


4.0 


0 


14 . 


4.47 


0^ 

0 

0 


15 . 


7.07 


8.14° 


16 . 


11,05 


5 . 2 ° 
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NUMERICAL RESULTS OF 16 POINT GRID 
COMPUTER SOLUTION 



POINT 


AMPLITUDE (cms.) 


PHASE (DEGREES 


1. 


3.47 


27.8 


2. 


3.43 


37.35 


3. 


3.16 


45.19 


4. 


3.16 


56.61 


5. 


4.60 


13.02 


6. 


4.94 


14.65 


7. 


4.95 


17.13 


8, 


3.74 


32.58 


9. 


• 

00 


49.20 


10. 


5.72 


9.11 


11. 


6.37 


8.11 


12. 


8.28 


3.06 


13. 


3.04 


2.98 


14. 


5.54 


5.41 


15. 


6.49 


7.76 


16. 


6.98 


7.30 
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